Longitudinal assessment of interstitial lung abnormalities on CT in patients with COPD using artificial intelligence-based segmentation: a prospective observational study

Background Interstitial lung abnormalities (ILAs) on CT may affect the clinical outcomes in patients with chronic obstructive pulmonary disease (COPD), but their quantification remains unestablished. This study examined whether artificial intelligence (AI)-based segmentation could be applied to identify ILAs using two COPD cohorts. Methods ILAs were diagnosed visually based on the Fleischner Society definition. Using an AI-based method, ground-glass opacities, reticulations, and honeycombing were segmented, and their volumes were summed to obtain the percentage ratio of interstitial lung disease-associated volume to total lung volume (ILDvol%). The optimal ILDvol% threshold for ILA detection was determined in cross-sectional data of the discovery and validation cohorts. The 5-year longitudinal changes in ILDvol% were calculated in discovery cohort patients who underwent baseline and follow-up CT scans. Results ILAs were found in 32 (14%) and 15 (10%) patients with COPD in the discovery (n = 234) and validation (n = 153) cohorts, respectively. ILDvol% was higher in patients with ILAs than in those without ILA in both cohorts. The optimal ILDvol% threshold in the discovery cohort was 1.203%, and good sensitivity and specificity (93.3% and 76.3%) were confirmed in the validation cohort. 124 patients took follow-up CT scan during 5 ± 1 years. 8 out of 124 patients (7%) developed ILAs. In a multivariable model, an increase in ILDvol% was associated with ILA development after adjusting for age, sex, BMI, and smoking exposure. Conclusion AI-based CT quantification of ILDvol% may be a reproducible method for identifying and monitoring ILAs in patients with COPD. Supplementary Information The online version contains supplementary material available at 10.1186/s12890-024-03002-z.


Background
Smoking and aging cause a variety of chronic respiratory diseases, such as chronic obstructive pulmonary disease (COPD) and interstitial lung disease (ILD).COPD is pathologically characterized by alveolar destruction, known as emphysema, and small airway disease, leading to persistent airflow limitations and respiratory symptoms [1], while ILD is characterized by inflammation and fibrotic changes in the lung parenchyma [2].Although the underlying pathophysiology of ILD is different from that of COPD, the coexistence of emphysema and ILD is common and known as combined pulmonary fibrosis and emphysema (CPFE) [3,4].This condition is associated with a high incidence of lung cancer and pulmonary hypertension, leading to a poor prognosis [3,5], but its pathogenesis remains unknown.
In contrast with histological analyses in which coexistence of emphysema and fibrosis is easily observed [6], fibrotic changes in emphysematous regions are subtle and cannot be resolved on high-resolution computed tomography (HRCT).However, interstitial lung abnormalities (ILAs) can be incidentally identified on CT in a subgroup of patients with COPD [7].The presence of ILAs in patients with COPD has a negative impact on clinical outcomes [8,9], and indeed, it is possible that ILAs are a precursor of CPFE.Therefore, the early detection of ILAs and assessment of their progression, particularly using an objective and reproducible quantification method, are crucial.
As the need for computer-aided CT analysis for ILD has increased, various methods of quantitative CT analysis have emerged, such as densitometry, texture analysis, and artificial intelligence (AI)-based methods [10].AI-based image analysis has advantages in terms of high reliability and reproducibility over other methods.Furthermore, indices derived from AI-based image analysis have been associated with pulmonary function decline, exacerbation, and prognosis in ILD [11][12][13].Artificial intelligence based quantitative CT image analysis software (AIQCT) is an AI-based image analysis software that can automatically classify chest CT images into normal lung, ground-glass opacities (GGOs), reticulations, consolidations, honeycombing, nodules, hyperlucencies, interlobular septum, bronchi, and vessels [14].Both normal lung regions and the bronchus volume quantified using AIQCT have been associated with a poor prognosis in patients with idiopathic pulmonary fibrosis (IPF) [14].However, the utility and applicability of AIQCT for quantifying ILAs in patients with COPD remains unclear.
Therefore, this study aimed to test the utility of AIQCT in evaluating the presence of ILAs in two independent COPD patient cohorts as well as ILA development using the longitudinal data from one of the cohorts.
Furthermore, the study explored abnormal CT patterns associated with the emergence of ILAs.

Study population and longitudinal analysis
This is a retrospective analysis of two prospective observational cohorts of patients with COPD [15][16][17][18].The discovery cohort included subjects who were older than 40 years, had a smoking exposure of more than 10 pack-years and were diagnosed with COPD according to the GOLD guidelines (forced expiratory volume in 1 s (FEV 1 )/forced vital capacity (FVC) less than 0.7 with respiratory symptoms) at Kyoto university hospital in 2006 to 2012.The exclusion criteria were (1) a diagnosis of ILD, (2) a prior history of lobectomy, and (3) inadequate CT images for quantitative analysis due to abnormal shadows such as pleural effusion, pneumonia, or poor quality.In addition to cross-sectional analysis, longitudinal analysis was performed on those who underwent chest CT scans within 5 ± 1 years from baseline (shown in Fig. 1).The validation of the thresholds of AIQCT indices for detecting ILA was performed using an independent dataset from another cohort [18].The validation cohort was a prospective observational cohort from Kyoto university hospital and a respiratory clinic (Terada Clinic), consisting of smokers aged ≥ 40 years with a history of ≥ 10 pack-years between April 2018 and April 2020.The details of the validation cohort are described in the additional file 1.

CT acquisition and visual assessment of ILAs
In the discovery cohort, all chest CT scans were acquired at full inspiration using an Aquilion 64 (sharp kernels, 0.5 mm slice thickness, 120 kVp, and autoexposure control; Cannon Medical Systems, Otawara, Japan).In the validation cohort, all chest CT scans were also acquired using an Aquilion Presicion scanner or an Aquilion lightning scanner (sharp kernels, 1.0 mm slice thickness, 120 kVp, and autoexposure control; Cannon Medical Systems, Otawara, Japan).Among one chest radiologist with 16 years' experience and two CT-experienced pulmonologists with ≥ 8 years' experience, two visually assessed ILAs in chest CT images according to the Fleischner Society statement [7].Any cases of discordance were settled by a chest radiologist or consultation between the two raters.The patients' clinical information was blinded, and the baseline and follow-up images were not compared when assessing the follow-up CT scans.

Quantitative CT analysis using AIQCT
Quantitative analysis of the chest CT scans was performed using AI-based image analysis software named AIQCT, which was originally developed in collaboration with FUJIFILM Corporation to quantify parenchymal and airway abnormalities in IPF [14].AIQCT uses a network architecture based on U-Net and automatically classifies structures on chest CT images into the following groups: normal lung, ground-glass opacities (GGOs), reticulations, consolidations, honeycombing, nodules, hyperlucencies, interlobular septum, bronchi, and vessels.According to a previous report [12], interstitial lung disease volume (ILDvol) was defined as the combined volume of honeycombing, reticulations, and GGOs.The absolute ILDvol and the percentage of ILDvol to total lung volume (ILDvol%) were analyzed in this study.

Statistical analysis
Data are presented as the mean ± standard deviation (SD) unless otherwise indicated.The kappa coefficient was calculated to assess interrater reliability in the visual assessment of ILA [19].Receiver operating characteristic (ROC) curve analyses and the Youden index were used to obtain the optimal threshold of ILDvol% in each lung volume fraction for detecting ILA.The areas under the ROC curves (AUCs) were compared using the DeLong test, and p values were corrected by the Bonferroni method.To test whether the diagnostic performance could be improved, sensitivity analysis was done in which the lung volume was equally divided into three or ten parts along the craniocaudal axis, and the highest ILDvol in each fraction was used in ROC curve analysis.Longitudinal changes in AIQCT indices were compared between subjects with new ILAs and those consistently without ILAs using the t test.Furthermore, adjusted mean differences in AIQCT indices were calculated, including age, sex, BMI, smoking status, and cumulative smoking exposure as covariates.Statistical analyses were performed using R statistical software version 4.3.1.

Study population and quantified CT features
Among 234 subjects, 32 (14%) had ILAs in the discovery cohort, and 15 out of 153 patients (9.8%) had ILAs in the validation cohort (Table 1).Interrater reliability for the visual assessment of ILAs was fair to moderate (kappa coefficient = 0.323-0.525,agreement = 81-88%).Compared to subjects without ILAs, those with ILAs were older and showed milder airflow limitation and a deteriorated diffusion capacity (Table 2).Table 2 also shows the AI-based quantification of CT features.The subjects with ILAs had higher incidences of honeycombing, reticulations, and GGOs than those without ILAs.Furthermore, nodule and bronchi volumes were greater in those with ILAs than those without ILAs.The clinical and radiological characteristics of the validation cohort are described in Table S1.Representative images of the AIQCT analysis are shown in Fig. 2.

Diagnostic performance of AIQCT indices
By using ROC curve analysis, the optimal threshold of ILDvol% on whole lung for identifying ILAs was determined.As shown in Figure S1, the optimal threshold of ILDvol% on whole lung was 1.203%, and the area under the ROC curve was 0.863, with 87.5% sensitivity and 77.2% specificity.In sensitivity analysis using ILDvol% within a third or tenth of the whole lung volume, the sensitivity and specificity were also high, but the diagnostic performance was not improved (shown in Fig. S1).When the optimal threshold of ILDvol% on whole lung was applied to the validation cohort, the diagnostic performance was good (sensitivity was 99.3% and specificity was 76.3%).

Longitudinal changes in the visual assessment of ILAs and AIQCT indices
124 patients were included in the longitudinal analysis who took follow-up CT scans in 5 ± 1 years after the baseline scans (Fig. 1).Eight (7%) patients demonstrated new ILAs on the follow-up CT (mean CT scan interval = 4.9 ± 0.5 years) (Fig. S2).When comparing the baseline characteristics of patients consistently without ILAs and those with new ILAs (Table 3), pulmonary function and AIQCT indices were not different between the groups.Figure 3 shows the longitudinal changes in the AIQCT indices.The ILDvol and the volumes of its components (honeycombing, reticulations, and GGOs) increased in patients with new ILAs at the follow-up CT scan.The volumes of nodules also increased in patients with new ILAs at the follow-up, which may represent an increase in the number of nodules.The volume changes in these indices were greater in patients with new ILAs than in those consistently without ILAs after adjusting for age, sex, BMI, and smoking exposure (Table 4).Bronchial volume tended to increase in patients with new ILAs, although the difference did not reach statistical significance (Table 4).Representative images illustrating the ILAs monitored during the follow-up period are shown in Fig. 4.
When using the optimal ILDvol% threshold defined in the cross-sectional analysis (1.203%), 6 out of 8 patients were classified as having no ILAs at baseline, but all were classified as having ILAs at follow-up.

Discussion
This study showed that AI-based quantitative analysis identified COPD with ILAs with high sensitivity and specificity.Furthermore, ILAs appeared in 7% of patients  with COPD 5 years after baseline CT acquisition.The volumes of not only GGOs, reticulations, and honeycombing but also nodules increased in patients with new ILAs.This is the first study to classify the findings from chest CT scans into detailed lesions, including airway lesions and to quantitatively analyze and longitudinally assess patients with COPD and ILAs.The prevalence of ILAs in patients with COPD in this study was 10 to 14%, which was consistent with the previous reports [20,21].ILAs in patients with COPD are more symptomatic, frequently exacerbated, and have a poor prognosis [8,20,22].Therefore, the effort to detect ILAs in patients with COPD is crucial for understanding and managing disease progression.Quantitative interstitial features assessed by local histogram-and distancebased method are associated with pulmonary functions, respiratory symptoms, and mortality in patients with COPD [23].The increase of quantitative interstitial features is associated with an annual decline in FVC, a decrease in exercise tolerance, and a poor mortality in longitudinal assessments [11].Furthermore, the progression of interstitial changes quantified using deep learning-based method had further impact on mortality in addition to the progression of emphysema in the COPD cohort [24].This study extends those findings by providing a longitudinal, quantitative assessment of patients with new ILAs via a detailed classification of chest CT scans using AIQCT.
AIQCT achieved good diagnostic performance with high sensitivity and specificity in identifying ILAs.This result is comparable to that from previous reports, although the quantifying techniques differed [25][26][27].The threshold values in this study (1.203%) and previous reports (1.8 to 3.6% [25][26][27] are also lower than the threshold value for visual assessment proposed by the Fleischner Society (5%) [7].The visual quantification of ILD volume based on tomographic images might be difficult, and the volume itself is usually overestimated [27].In the sensitivity analysis, the diagnostic performance was not improved even when the lung volume was divided.This might be because ILD lesions are usually localized in the basal part of the lung in patients with ILAs, and artifacts on CT images did not affect the diagnosis of quantitative ILAs.
ILAs emerged in 7% of patients with COPD within 5 years in this study.The prevalence of new ILAs was higher than that in a population-based study [28].This study included more at-risk patients, that is, older patients and those with greater smoke exposure [29].The baseline pulmonary function and volumes of CT regions did not differ between patients who developed ILAs and those who consistently remained without ILAs, which indicates that the occurrence of ILAs cannot be predicted by baseline pulmonary function or AIQCT indices.
The longitudinal changes in volumes of ILD regions and nodules were greater and that in bronchial volumes tended to be greater in patients with new ILAs than in those consistently without ILAs.Recently, some studies have shown that small airways are important to the pathogenesis of ILD or ILAs [30][31][32][33].The microCT analysis of unbiased, systematic uniform random tissue samples from patients with IPF showed that the terminal bronchiole was dilated even in the presence of lesions without fibrosis [31].The airway wall on HRCT is thicker in subjects with ILAs than in those without ILAs [30].Furthermore, bronchiectasis and peribronchial fibrosis were identified in ILA lungs in histopathological studies [32].These airway changes may appear as nodules due to the resolution of HRCT.This study supports the hypothesis that interstitial fibrotic changes may start in the bronchial regions.
The interrater agreement was not high, neither between CT-experienced pulmonologists or with respect to a chest radiologist.Airspace enlargement with fibrosis, scar-like lesions, and paraseptal emphysema may make it harder to assess the existence of ILAs, especially in patients with COPD [34].AIQCT is useful for providing objective and reproducible assessments of even slight interstitial abnormalities in emphysematous patients.Furthermore, it can detect the progression of ILAs.To elucidate the pathogenesis of ILAs and their clinical sequences to ILD, repeated histological analyses are  Data are the mean ± standard deviation or adjusted mean differences with 95% confidence intervals in parentheses.The differences between patients who did not develop ILA and those who developed ILA were adjusted by age, sex, BMI, smoking status, and cumulative smoking exposure.ILA interstitial lung abnormality; ILD interstitial lung disease needed.However, biopsy is too invasive to repeat over time, and the biopsy procedure may alter the microenvironment of the lung.In contrast, CT image analysis is less invasive and allows repeated examinations without changing the lung microenvironment.Imaging studies using AIQCT may elucidate the pathogenesis of ILAs.This study has some limitations.First, the number of new ILAs was small, making it is difficult to explore factors associated with the future development of ILAs.Second, a substantial number of patients did not undergo follow-up CT scans within the period.A selection bias may exist when patients undergo follow-up CT scans.Third, most patients included in this study were male because COPD patients tend to be male in Japan and Asian countries [35,36].

Conclusions
In conclusion, this study showed that ILAs in patients with COPD could be measured using AIQCT.7% of patients with COPD showed new ILAs at follow-up, and these patients could be identified by using AIQCT.Although additional studies in longitudinal cohorts are needed, AIQCT revealed that honeycombing, reticulation, GGO, and nodule volumes were greater in patients who developed ILAs.

Fig. 2
Fig. 2 CT images of COPD with ILAs and corresponding AIQCT images Legend: The axial CT image of a 68-year-old male patient with COPD and ILAs in the discovery cohort (Fig. 2. a) and the corresponding AIQCT image overlaid on the original image (Fig. 2. b).Sagittal image of the same patient (Fig. 2. c) and the corresponding AIQCT image (Fig. 2. d).The ILD volume was 150 mL (2.90%)

Table 1
Patient characteristics in the two cohorts Data are the mean ± standard deviation or the number of patients with percentage in parentheses.BMI body mass index; mMRC> modified British Medical Research Council Dyspnea scale; FVC forced vital capacity; FEV 1 , forced expiratory volume in 1 s; ILA interstitial lung abnormality

Table 2
Clinical and radiological comparisons between patients with and without ILA in the discovery cohort 1, forced expiratory volume in 1 s; DLCO diffusing capacity of the lung for carbon monoxide; ILD interstitial lung disease

Table 3
Baseline characteristics of patients with COPD without ILAs at enrollment who subsequently developed ILAs and those who consistently showed no ILA findings on follow-up CT scans Data are the mean ± standard deviation or the number of patients with the percentage in parentheses.ILA interstitial lung abnormality; BMI body mass index; mMRC modified British Medical Research Council Dyspnea scale; FVC forced vital capacity; FEV 1 , forced expiratory volume in 1 s; DLCO diffusing capacity of the lung for carbon monoxide; ILD interstitial lung disease

Table 4
Longitudinal changes in AIQCT indices